Conformal cloak for waves 
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Conformal invisibility devices are only supposed to work within the validity range of geometrical 
optics. Here we show by numerical simulations and analytical arguments that for certain quantized 
frequencies they are nearly perfect even in a regime that clearly violates geometrical optics. The 
quantization condition follows from the analogy between the Helmholtz equation and the stationary 
Schrodinger equation. 
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Invisibility came into sight as the first nontrivial appli- 
cation of transformation optics [T]-[5] . The progress made 
has been impressive, but complete cloaking devices have 
never been demonstrated in practice yet. All electromag- 
netic cloaking experiments reported in the literature were 
for reduced cases. For example, the first cloaking device 
[6] worked for microwaves of one frequency and polariza- 
tion, and this only in the approximation of geometrical 
optics. Cloaking in tapered waveguides is also approx- 
imative [7J. Carpet cloaking [8, 9 is a drastic form of 
reduced cloaking where an object is not made to disap- 
pear completely, but to appear as being flat. There are 
two reasons for the need of resorting to reduced cloaking 
devices, a practical and a fundamental one. The ma- 
terials required for perfect cloaking [ll |3] are extremely 
difficult to fabricate, because they need to implement 
impedance-matched anisotropic media [5 , and therefore 
perfect cloaking is impractical. The fundamental prob- 
lem is that perfect cloaking also implies propagation with 
a superluminal phase velocity that reaches infinity [10], 
which is possible in principle, but only for discrete fre- 
quencies. Perfect cloaking in a broad band of the spec- 
trum is therefore physically impossible. Note that these 
qualifications do not apply to acoustic cloaking [TTl [12] 
where the required materials are much easier to manu- 
facture and the constraints from relativistic causality are 
not relevant. 

One of the earliest ideas for invisibility devices. Optical 
Conformal Mapping [2 , 14 , has the advantage of requir- 
ing rather ordinary optical materials — they are opti- 
cally isotropic and nonmagnetic (but still need a large 
index range). It can be easily applied in acoustics, be- 
cause it does not rely on materials with artificially de- 
signed anisotropic mass [TTl IE] • This cloaking method 
was originally derived [2] for the approximation of ge- 
ometrical optics (or acoustics). In isotropic media, per- 
fect invisibility is mathematically impossible, not only for 
physical reasons, because the inverse scattering problem 
is uniquely solvable [13] and so waves cannot completely 
hide from the fact that they propagate in media. Here 



we show, however, that conformal cloaking can be nearly 
perfect for discrete frequencies. 

Cloaking by Optical Conformal Mapping consists of 
two parts: (1) the implementation of a conformal coordi- 
nate transformation and (2) an index profile with certain 
properties that cannot be reduced to flat space by a co- 
ordinate transformation, a non-Euclidian part [TSl [16]. 
A conformal coordinate transformation is valid for com- 
plete waves, but it cannot make a cloaking device on its 
own [2'; the non-Euclidean part completes the device, 
but it is only supposed to work within the validity range 
of geometrical optics. Full-wave simulations of Optical 
Conformal Mapping are rare; in a recent review [17 the 
effect of conformal coordinate transformations was stud- 
ied, but without the non-Euclidean part. Not surpris- 
ingly, cloaking is not possible here, but it is wrong to 
conclude from an incomplete numerical experiment that 
cloaking by Optical Conformal Mapping would not work. 
The simulations we show here prove to the contrary, even 
in a regime far away from geometrical optics. 

Let us begin with a brief recapitulation of Optical Con- 
formal Mapping [2^, . Consider a planar medium with 
the two-dimensional graded index profile n{x^y). The 
medium is purely electrical such that e = n'^. In this 
case, electromagnetic waves polarized such that the elec- 
tric field points in vertical direction (orthogonal to the 
plane) obey the Helmholtz equation [8] 

^={dl^dy^ n^k^) ^ = (49.*^. + ^'^') ^ (1) 

where i/j denotes the electric-field component, z = x -\-iy 
and k is the wavenumber. Electromagnetic waves polar- 
ized such that the magnetic field points in vertical di- 
rection also obey the Helmholtz equation Q, but only 
approximately in the regime of geometrical optics [5 . In 
this case denotes the magnetic- field component. A lin- 
ear combination of the two polarizations constitutes an 
arbitrary electromagnetic wave in the planar medium. 
Under a conformal mapping w = w{z) the Helmholtz 
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equation in w space appears as 

{4d^d^^n'^k^)ilj = 0, n = n' 



dw 



dz 



(2) 



In the language of transformation optics [U |5] iL'-space 
is called virtual space and z-space is the physical space. 
Consider the simplest mapping suitable for cloaking, the 
Zhukowski transform O [5] 



z 



or, equivalent ly, z 
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that maps a virtual i(;-space with two Riemann sheets 
on two regions of physical z-space, one is the exterior 
and the other the interior of a circle of radius a (Fig. 1). 
If the virtual space were empty light rays would travel 
along straight lines in K;-space that may cross the branch 
cut from the exterior to the interior w sheet and then get 
absorbed at the infinity of the latter that corresponds to 
the singularity of w{z)^ which, for the Zhukowski trans- 
form (|3|, lies at z = in physical space. This feature has 
caused the distortion of wave propagation [17 mentioned 
above. However, it is possible to shepherd the lost rays 
back to the exterior sheet by filling the interior sheet with 
a suitable index profile that cannot be reduced to empty 
space by a coordinate transformation [2l [14]. Many of 
such profiles are possible [5] , including ones with negative 
refraction ^ (that do not cause phase delays). What 
these index profiles have in common is their ability to 
let light propagate in closed trajectories. Consider the 
simplest case, the Hooke and the Kepler profile [2| [14]: 
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(4) 
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where Wi is the branch point 2a of the Zhukowski map 
(|3| and ro is chosen to be 4a. The index n in physi- 
cal space ranges from to 14.5 for the cloaking device 
with Hooke profile and from to 13.3 for the Kepler 
case. In both cases, the ray trajectories are closed curves 
in virtual space (in particular ellipses) such that after a 
loop in the "underworld" — on the interior sheet — the 
light returns to the exterior sheet and propagates along 
straight lines in iL'-space in the same direction it came 
from. In the "underworld" the light is confined to a circle 
around wi with radius ro, because beyond this circle n' 
is purely imaginary and so light cannot propagate there 
(light waves decay exponentially). In physical space, the 
exterior of the virtual circle is the interior of its map z{w)^ 

i.e. of the curve r = a ^2 + cos (j) — y^(2 + cos^)^ ~ ) 
polar coordinates. This region — and anything inside it 
— is invisible, while light bends around it such that rays 
are asymptotically straight lines in physical space: the 
conformal transformation ([3| combined with one of the 




FIG. 1: In transformation optics, electromagnetic waves are 
transformed from an empty virtual space to physical space 
by an appropriate medium. In Optical Conformal Mapping 
this is achieved by a conformal transformation that requires 
an optically isotropic material. The virtual space consists of 
Riemann sheets (above) that are mapped to the physical plane 
(below). However, light may cross the branch cut between two 
sheets and get absorbed at a singularity of the transformation, 
casting a shadow with zero amplitude, as the figure illustrates 
for the Zhukowski map (|3]). Cloaking with optically isotropic 
materials (Fig. 2) is possible when the lower sheet contains a 
medium where light propagates in closed loops (Fig. 3). 

profiles Q and ([5| makes a cloaking device. Note that 
light rays are refracted at the interface between the in- 
terior and the exterior region of the Zhukowski map (|3| 
where the profiles Q and ([5| reside. Rays are refracted 
back to their original direction after one loop on the in- 
terior sheet, but waves are partially refiected at sudden 
index steps. When the index varies much more rapidly 
over the scale of the wavelength, as it is the case for 
index steps, the conditions of geometrical optics [5 are 
violated, causing the conversion of one wave front into 
two, i.e. partial refiection. Furthermore, the detour in 
the "underworld" causes a uniform time delay [14 that 
appears as a phase shift in the part of the wave that en- 
tered there, and thus creates phase dislocations that may 
cause diffraction overshadowing the cloaking of small de- 
vices. 

Here we have tested the performance of conformal 
cloaking devices for waves using simulations made with 
standard commercial software. In the invisible region 
of the device we used an imaginary index profile. One 
can also shield this region with perfect mirrors; our sim- 
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ulations show little principal difference between the two 
(only quantitative differences). We also introduced a cut- 
off radius Tc in physical space that makes our cloaking de- 
vice finitely extended. Conformal transformations such 
as the Zhukowski map ([3| act across the entire plane, 
in contrast to quasiconformal transformations 0, so the 
medium implementing the transformation is, in princi- 
ple, infinitely extended. Yet beyond the cut-off radius 
Vc = 5a we put n = 1 and observed that this cut-off 
has a negligible effect on the propagation of waves. More 
importantly, our simulations show that the sharp index 
steps at the interface of the inner cloaking region do nei- 
ther cause reflections nor phase delays for certain discrete 
wavenumbers of light (Fig. 2). 




FIG. 2: Color online. Simulations of wave propagation for a 
cloaking device based on the Zhukowski map (Is]) and the Ke- 
pler profile ([5]); our results for the Hooke pro file (|4| are very 
similar. In each panel the outer circle describes the boundary 
of the device (with cut-off radius rc), the inner circle contains 
the core of the cloaking device that carries the transformed 
profile ([5]) and the pupil indicates the cloaked region where 
the wave decays exponentially. The top pictures (a,b) show 
the cloaking of incident plane waves and the bottom pictures 
(b,d) cylindrical waves. Our simulations illustrate two ex- 
treme cases, kro = 5 (a,c) with nearly perfect invisibility and 
kro — 6 (b,d) with a pronounced phase dislocation and a re- 
sulting diffraction pattern. Note that the core of the device is 
small in comparison with the wavelength such that geomet- 
rical optics is no longer a good approximation, yet cloaking 
can be nearly perfect. 



The wavenumbers of nearly-perfect invisibility turn 
out to be related to the eigen-frequencies of light in the 
index profiles Q and ([5|. We can easily deduce them 
from the analogy between the Helmholtz equation Q in 
virtual space and the stationary Schrodinger equation 



2m , 



(^-/7) U = 0, E-U 



2m 



(6) 



where U denotes the potential, E the energy and m 
the mass. The Hooke profile Q corresponds to the 2D 
harmonic-oscillator potential 



— 2" ^ ~ ^1 



(7) 



with ujo being the oscillation frequency. An eigenstate of 
the 2D oscillator has the quantized energy 
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for non- negative integer I. From relations ^ and ([8| 
follows 



kro = 2(/ + 1) (Hooke) 



(9) 



that defines the wavenumber k of an eigenmode. Along 
similar lines we deduce the eigenmodes in the Kepler pro- 
file ([5| from the 2D hydrogen spectrum [19 , and obtain 
in this case 



kro = 2/ + 1 (Kepler). 



(10) 



It is clear from semiclassical quantum mechanics that for 
these eigenmodes the phase difference along the closed 
loop on the interior sheet in virtual space is an integer 
multiple of 27r. Therefore the device does not cause a 
phase dislocation; the emerging wavefront on the exte- 
rior sheet is intact. What seems more surprising is the 
fact that the incident wave is not reflected at the bound- 
ary between exterior and interior sheet (within the ac- 
curacy of our simulations) even when the index jumps 
by an infinite amount, as is the case for the Kepler pro- 
file ([5| in virtual space. Finally, we observed (Fig. 2) 
that the cloaking device performs nearly perfectly for the 



wavenumbers and (10), showing no signs of reflection 



and diffraction beyond free wave propagation. Cloaking 
is almost perfect even when the device is smaller than the 
wavelength of light, a regime that lies outside the validity 
range of geometrical optics [5 . 

The perfect transmission and hence perfect invisibility 
at an eigenfrequency can be made plausible by a simple 
one-dimensional model. Imagine, instead of the lower 
Riemann sheet (Fig. 1) where the light enters and leaves 
through the branch cut, a one-dimensional index profile 
(Fig. 3) that light enters from one side and leaves at 
the other. This simple model contains the essence of 
the more complicated behavior of light propagation on 
the interior Riemann sheet where the index jumps, but 
after one loop returns to the original value; in the model 
we have opened this loop such that the two sides of the 
branch cut are represented by two interfaces. We denote 
the reflectivities and transmittivities of each interface by 
ri,r2 and ^1,^2, and obtain from the Fresnel coefficients 
[5] the relations 



ri 



-^2 : 



tlt2 = 1 



1 



(11) 
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formations [3^ — and this even in a regime far beyond 
geometrical optics. 
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FIG. 3: Color online. One-dimensional model. The top pic- 
ture shows a typical closed loop of a light ray on the Riemann 
surface of the Zhukowski map ^ with the Kepler potential 
([5| on the lower sheet. The bottom picture visualizes the 
refractive index profile n experienced by the ray on its tra- 
jectory in virtual space. At the sharp interfaces between the 
medium a wave is reflected and transmitted with coefficients 
Vi and ti. However, when the phase delay cp in the dielectric 
structure is an integer multiple of 27v no net reflection occurs 
due to interferences at multiple reflections. In this case the 
structure is perfectly invisible. 

Consider a wave i/j propagating in x direction in a uni- 
form background medium with index no towards the non- 
uniform profile representing the interior sheet. Part of 
the wave is reflected, the transmitted part propagates 
to the second interface with phase (f and is multiply re- 
flected inside the index profile, such that in total 

^ ^ e^^°^^+e-^^°^^(ri+tie^^r2e^^pt2) , (12) 



^=0 



For an eigenmode the phase delay cp is an integer mul- 



tiple of 27r. In this case the Fresnel relations (11) imply 
that i/j = exp(inoA:x); the interference of light by multiple 
reflections causes no reflection at all — the wave is per- 
fectly transmitted and so the index profile is invisible. It 
is remarkable that this behavior carries over to the wave 
propagation on virtual Riemann sheets that, for certain 
frequencies, are not only invisible themselves but make 
everything inside invisible. 

We thus showed by numerical simulations and sim- 
ple analytical arguments that cloaking devices based on 
Optical Conformal Mapping [2 work perfectly for dis- 
crete frequencies — like the cloaking devices made by 
the implementation of non-conformal coordinate trans- 



[1] A. Greenleaf, M. Lassas, and G. Uhlmann, Math. Res. 

Lett. 10, 685 (2003); Physiol. Meas 24, 413 (2003). 
[2] U. Leonhardt, Science 312, 1777 (2006). 
[3] J.B. Pendry D. Schurig, and D.R. Smith, Science 312, 

1780 (2006). 

[4] H. Chen, C.T. Chan, and P. Sheng, Nature Materials 9, 
387 (2010). 

[5] U. Leonhardt and T.G. Philbin, Geometry and Light: 
The Science of Invisibility (Dover, Mineola, 2010). 

[6] D. Schurig, J.J. Mock, B.J. Justice, S.A. Cummer, J.B. 
Pendry A.F. Starr, and D.R. Smith, Science 314, 977 
(2006). 

[7] I.I. Smolyaninov, V.N. Smolyaninova, A.V. Kildishev, 

and V.M. Shalaev, Phys. Rev. Lett. 102, 213901 (2009); 

S. Tretyakov, P. Alitalo, O. Luukkonen, and C. Simovski 

Phys. Rev. Lett. 103, 103905 (2009). 
[8] J.S. Li and J. B. Pendry Phys. Rev. Lett. 101, 203901 

(2008). 

[9] R. Liu, C. Ji, J.J. Mock, J.Y. Chin, T.J. Cui, and D.R. 
Smith, Science 323, 366 (2009); J. Valentine, J.S. Li, 
T. Zentgraf, G. Bartal, and X. Zhang, Nature Materials 
8, 568 (2009); L.H. Gabrielh, J. Cardenas, C.B. Poitras, 
and M. Lipson, Nature Photonics 3, 461 (2009); T. Er- 
gin, N. Stenger, P. Brenner, J.B. Pendry, and M. We- 
gener, Science 328, 337 ( 2010); B. Zhang , Y. Luo, X. 
Liu, and G. Barbastathis, |arXiv:1012.2238| X. Chen, Y. 
Luo, J. Zhang, K. Jiang, J.B. Pendry, and S. Zhang, 
larXiv:10T2 '.2783 
[10] U. Leonhardt and T.G. Philbin, New. J. Phys. 8, 247 

(2006) . 

[11] H. Chen and C. T. Chan, Appl. Phys. Lett. 91, 183518 

(2007) . 

[12] S. Zhang, C. Xia, and N. Fang, Phys. Rev. Lett. 106, 

024301 (2011). 
[13] A.I. Nachman, Ann. Math. 128, 531 (1988). 
[14] U. Leonhardt, New. J. Phys. 8, 118 (2006). 
[15] U. Leonhardt and T. Tyc, Science 323, 110 (2009). 
[16] T. Tyc, H. Chen, C. T. Chan, and U. Leonhardt, IEEE 

Journal of Selected Topics in Quantum Electronics 16, 

418 (2010). 

[17] Y.A. Urzhumov, N.B. Kundtz, D.R. Smith, and J.B. 

Pendry J. Opt. 13, 024002 (2011). 
[18] T. Ochiai, U. Leonhardt, and J.C. Nacher, J. Math. Phys. 

49, 032903 (2008). 
[19] B. Zaslow and M.E. Zandler, Am. J. Phys. 35, 1118 

(1967). 



